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MECHANISM AND KINETICS OF CORROSION OF SELECTED IRON AND 
COBALT ALLOYS IN REFLUXING MERCURY 


by Louis Rosenblum, Coulson Scheuermann, Charles A. Barrett, 

and Warren H. Lowdermilk 

Lewis Research Center 
SUMMARY 

Five alloys of engineering interest, the cobalt -base HS-25 and H-8187 and the iron- 
base SICROMO-9M, AM-350, and AM-355, were corrosion tested in boiling mercury in 
the temperature range of 1000° to 1300° F (811° to 977° K) for times as long as 5000 
hours. Metallographic, chemical, and physical analyses were used to determine the ex- 
tent and the nature of the corrosion. 

Corrosion mechanism and reaction kinetics were inferred through a comparison of 
experimental and theoretical kinetic constants. The experimental kinetic constants (i. e. , 
the time-law constant b p and the activation energy AH r ) were obtained from the test 
data by multiple regression analysis. The theoretical kinetic constants b. and AH^ 
for each of i steps of proposed corrosion models were obtained from transport theory 
relations and the pertinent boundary conditions. 

Changes of mechanism during the progress of corrosion were indicated for all the 
alloys tested except the AM-350 and AM-355 alloys. The shifts from one corrosion re- 
gime to another were associated with changes in the nature and length of the corrosion- 
affected zone. Three corrosion regimes were characterized, linear I, parabolic, and 
linear II. The probable rate -determining step in the two linear rate regimes is 
boundary-layer diffusion. The probable rate -determining step in the parabolic rate re- 
gime is liquid diffusion. 


INTRODUCTION 


The present interest in space power systems that employ mercury as the thermo- 
dynamic working fluid at temperatures in the range of 1000° to 1200° F (811° to 922° K) 
(refs. 1 and 2) has prompted the need for an understanding of the mechanism and kinetics 
of mercury corrosion. Such knowledge could be an invaluable guide in establishing oper- 



ational limits, extrapolating performance, and translating small-scale test results to 
larger system applications. 

According to the prevailing concepts of reaction mechanisms, a reaction such as the 
corrosion of a solid material by mercury occurs as a sequence of steps. The combined 
effect of these steps determines the net rate of the reaction. However, if, as is common, 
one step is significantly slower than the others, the rate of the entire reaction is gov- 
erned by the slowest step. In cases where the slow step involves a diffusion process 
(e. g. , volume diffusion in a solid), the reaction is termed diffusion controlled and is 
governed by the laws of diffusion kinetics. Alternately, if a boundary-region process 
(e. g. , solution) constitutes the slow step, the rate of the overall reaction is determined 
by the kinetics of this process. 

The kinetics of each step of a corrosion mass -transport reaction can be expressed 
by the empirical relation 



(1) 


where w. is a measure of the mass transported, b- is the time-law constant for the 
th 1 1 

i step, and t is time. The quantity k- is known as the specific rate constant for the 
i^ step and is a function of temperature only; that is, 


k ‘ = a ‘ exp (i?) <2) 

where a. is a constant, R is the gas constant, T is absolute temperature, and AH. is 

1 th 1 

the ’’energy of activation’’ of the i step. 

Therefore, when a corrosion reaction is governed by one controlling step, the mass 
transport and the specific rate constant for the entire reaction will be equal to, respec- 
tively, the mass transport w g and the specific rate constant k g of the rate -controlling 
step. This relation, then, forms a basis for determining the mechanism and kinetics of 
corrosion. 

Only two mechanism studies of high-temperature dynamic mercury corrosion have 
been reported (refs. 3 and 4). In each study the same data were used, namely, the data 
for the corrosion of low -carbon steel in mercury over the temperature range of 900° to 
1200° F (755° to 922° K) presented in reference 5. Nevertheless, each study reached a 
different conclusion as to the mechanism of attack. The basic approach followed can be 
summarized as follows: 

(1) An assumption was made of the rate -controlling step. 

(2) A heat-transfer - mass-transfer analog equation was derived for k . 
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(3) Estimates were made of the mass-transfer coefficients and empirical constants 
of the heat-transfer - mass-transfer analog equation and a solution was obtained for k g . 

(4) The experimentally observed corrosion rate was compared with k for agree- 

s 

ment. 

(5) The activation energy derived from the experimental data was compared for 
agreement with the activation energy estimated to be associated with the assumed rate- 
controlling step. 

In the first study (ref. 3), liquid-boundary-layer diffusion was assumed to be the 
rate-controlling process. Order -of -magnitude agreement was obtained between the rate 
predicted by the derived heat-transfer - mass-transfer analog equation and the observed 
corrosion rate. The order-of -magnitude agreement between observation and theory was 
considered to justify the general validity of the premises made. 

A second study (ref. 4) used the same data and mass-transfer equation as the first 
study. It was concluded therein that the difference between the estimated and experimen- 
tal activation energies and the disparity between the observed and theoretical corrosion 
rate suggested that a process other than diffusion- control solution attack (probably inter- 
granular attack) is the controlling factor in the corrosion of steel by mercury. 

For completeness, corroboration of a corrosion mechanism by this approach re- 
quires the formulation of a number of corrosion-rate equations equal to the total number 
of possible rate -determining steps. Concomitantly, values of the pertinent material- 
transport coefficients and empirical constants are required to solve each equation. In 
many instances, these values have not been determined or are not known accurately. 

In the present study, a somewhat different approach was used. No assumption of a 
rate -controlling step was made. Rather, rate and rate -controlling step were inferred 
from a comparison of experimental and theoretical kinetic constants. The major prob- 
lem of this approach is the paucity of available data needed in the estimate of activation 
energy for several pertinent processes. In brief, the approach was as follows: 

(1) Values of the experimental kinetic constants b r and AH r were obtained from a 
multiple linear regression fit of the corrosion data to the estimating equation 

AH 

log p = a + b log t + (3) 

r r r 2.303 RT 

where p r is observed maximum corrosion penetration or wall recession (assumed to be 
proportional to w r ), t is test time, T is test temperature, and a r is a constant. 
Equation (3) was derived by the elimination of k between equations (1) and (2) (where the 
subscript r replaces i). 

(2) A corrosion model was constructed consistent with the general nature of the ob- 
served corrosion attack. For each step of the corrosion model, an equation was derived 
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(from considerations of transport theory and the relevant boundary conditions) in the 
form of equation (1), which yielded bj directly. 

(3) Values of AH^ were obtained from the relation 


AH.= 


R d(ln kj) 


d Ik 


( 4 ) 


which is the differential form of equation (2). 

(4) The theoretical kinetic constants AH^ and the b^ for each step in the proposed 
corrosion mechanism were compared with AH and b„, respectively. Agreement of 

i. r * 

the theoretical kinetic constants for the i in step with the respective experimental kinetic 

ii. 

constants was taken to infer that the i in step is rate determining. 

The method of data analysis employed in this approach provides a ready opportunity 
to test for the possibility of mechanism changes during the process of corrosion. Such 
changes may be expected, inasmuch as corrosion induces physical changes in the test 
material which can lead to changes in the material-transport path and, thereby, in the 
rate-controlling step. This process is analogous to certain oxidation reactions where 
changes in oxidation kinetics accompany changes in the thickness or adherence of oxide 
layers (ref. 6). The possibility of mechanism changes during the course of corrosion 
was explored. First, the data for each test material were grouped on the basis of the 
physical appearance of the corroded area. Second, a multiple linear regression fit of the 
grouped data to equation (3) was made. Third, a statistical goodness -of -fit comparison 
was made between the grouped and nongrouped data. 

For the experimental part of this study, 108 reflux capsule tests were performed on 
five materials of engineering interest, AM-350, AM-355, SICROMO-9M, HS-25, and 
H-8187. The test temperature range was 1000° to 1300° F (811° to 977° K) and test 
times were as long as 5000 hours. Metallographic analyses were performed on all test 
specimens to determine corrosion penetration. In addition, chemical and physical analy- 
ses were performed on selected specimens to determine the composition of the corrosion 
zones, the corrosion deposits, and the bulk alloys. (The authors express their apprecia- 
tion to Dr. Leo F. Epstein of General Electric Vallecitos Laboratory and Dr. Malcolm 
Parkman of Varian Associates, Palo Alto, California for helpful review and for pointing 
out an error in activation -energy estimation in the original manuscript. ) 
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EXPERIMENT 


Procedure 

The effect of mercury corrosion on HS-25, H-8187, SICROMO-9M, AM-350, and 
AM-355 was determined by means of reflux capsule tests wherein the capsule wall served 
as the test specimen. Table I lists the composition of the test alloys and table n lists the 
forms in which they were received, their prior history, and their sources. 


TABLE I. - COMPOSITION OF TEST ALLOYS 


Alloy 

Ni 

Co 

Fe 

Cr 

Mn 

Mo 

W 

Other 






wt. % 



HS-25 

9.9 

51. 1 

2.5 

19.8 

1.2 

— 

15. 5 


H-8187 

1.2 

64.6 

1 . 1 

20.6 

.6 

0 . 1 

11.8 


SIC HOMO -9M 

— 

— 

90.3 

8.6 

.5 

.6 

— 


AM-350 

4. 1 

— 

75. 1 

16.6 

1 . 1 

2.9 

— 

0.2 (N+C) 

AM-355 

4. 3 

— 

76.2 

15.5 

1.0 

00 

c\i 

— 

.2 (N + C) 


TABLE II. - FORM, HISTORY, AND SOURCE OF TEST ALLOYS 


Alloy 

Form 

History 

Source 

HS-25 

1/2-in. (1. 27-cm) round 
stock, centerless ground 

Solution treated at 
2250° F (1505° K) 

Haynes Stellite Co. 

H-8187 

2-in. (5. 08-cm) hot- 
worked square 

Solution treated at 
2250° F (1505° K) 

Haynes Stellite Co. 

SIC HOMO -9M 

1-in. (2. 54-cm) round 
stock 

Annealed at 1650° F 
(1171° K), cooled 
50° F (28° K) 
per hour 

Timken Roller 
Bearing Co. 

AM-350 

9/16-in. (1. 43-cm) rod, 
centerless ground 

Annealed at 1425° F 
(1047° K), air cooled 

A lie gheny - Ludlum 
Steel Corp. 

AM-355 

9/16-in. (1. 43-cm) rod, 
centerless ground 

Annealed at 1025° F 
(825° K), air cooled 

Allegheny - Ludlum 
Steel Corp. 


5 



Figure 1. - Reflux capsule. 


Capsule preparation. - Capsules (fig. 1) were machined from the as-received stock 
to the nominal dimensions of 1/2 -inch (1. 27-cm) outside diameter, 1^ inches (4. 45 cm) 
in length and 0. 040-inch (0. 102-cm) wall thickness. The inner surfaces were lapped to 
a mirror-like finish. The capsules were then degreased and chemically cleaned. 

Next, the capsules were filled to approximately one-third their volume with film- 
free, triple -distilled mercury and sealed by electron-beam welding in a vacuum of about 
10”° torr (10”° N/m' 5 ). During the welding process, which took less than 8 seconds, the 
capsules were held in a water-cooled jacket to prevent vaporization of the mercury. 

Heat treatment of the capsule as a result of the weld was confined to the vicinity of the 
weld. The sealed capsules were inspected radiographically to check the mercury level 
and the degree of wetting of the capsule wall. 

Capsule testing. - The capsules were tested in argon. Heating of the capsules was 
accomplished by means of individual capsule furnaces which supplied heat to the lower 
two-thirds of the capsule. The furnaces consisted of a Nichrome wire winding on a mag- 
nesium silicate core surrounded on the sides and bottom with 0. 375 inch (0. 952 cm) of 
magnesium silicate. The tops of the capsules were exposed so that heat was lost to the 
environment by natural convection and radiation. 

Test beds of 10 furnaces were employed (fig. 2). Each furnace was individually 
controlled to within ±10° F (±5. 5° K) of a preset temperature by its own powerstat. 
Chromel-Alumel thermocouples spot welded to the top and bottom capsule surface were 




Figure 2. - Test bed used for reflux capsule testing. 


1 C-55896 


used to measure temperature. The thermocouple output was continuously recorded on a 
multipoint strip -chart recorder. 

Tests were conducted at 1000°, 1100°, 1150°, 1200°, 1250°, and 1300° F (811°, 

866°, 894°, 922°, 950°, and 977° K); test times were as long as 5000 hours. The tem- 
perature indicated by the bottom thermocouple was used to represent the test tempera- 
ture. 

Post-test processing . - At the conclusion of the tests, the capsules were allowed to 
cool to room temperature, pierced by drilling, and drained of the mercury. The cap- 
sules were cut lengthwise along the axis, and residual mercury was removed by vacuum 
evaporation. 

One section of each capsule was mounted and examined metallographically in the 
etched and unetched conditions. Depths of corrosion penetration were measured with a 
filar eyepiece. Specimens of particular interest were analyzed by means of an electron- 
beam probe microanalyzer to determine the elemental composition of the corroded areas. 
Selected corrosion deposits were analyzed chemically to determine their composition. 


Results 

Test data . - A summary of test conditions, maximum corrosion-penetration depths, 
and appearance of the corrosion-zone structure is given in table m. 
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TABLE m. - REFLUX CAPSULE TEST DATA 


Capsule 

Temper- 

ature 


Test 

time, 

hr 

Maximum depth of attack 

Corrosion manifested as - 


Capsule 

Temper- 

ature 


Test 

time, 

hr 

Maximum depth of attack 

Corrosion manifested as - 

in. 

m 

in. 

m 

°i 

F 

°K 

°F 

°K 

SIC ROMO -9M 








AM 

-355 



1 

1000 

811 

308 

0.05X10' 3 

0. 13X10' 5 

Grain-boundary penetration 


1 

1100 

866 

1004 

0.34X10" 3 

0. 86X1CT 5 

Wall recession 

2 





308 

. 12 

.31 



2 

1100 

866 

1005 

.44 

1.12 



3 





1007 

.49 

1.24 



3 

1200 

922 

1004 

.83 

2.11 



4 





1007 

.51 

1.30 



4 

1200 

922 

1004 

.68 

1.73 



5 





2022 

.39 

.99 



5 

1200 

922 

1005 

.90 

2.29 



6 





4993 

1.28 

3.25 









o 

c 



7 





4993 

1.72 

4.37 


r 

6 

1300 

977 

1004 

1. 12X10"' 3 

2.84X10" & 

Wall recession 












in 



i nri/i 

1. 17 










** 



( 

lol. 

JU 

y < 

f 


z. y i 

wan re 

cession 

8 

1100 

866 

305 

0. 17X10 J 

0.43X10 0 

Grain -boundary penetration 











9 





305 

.29 

.73 












10 





1005 

.90 

2.29 



1 

1100 

866 

303 

1.29X10' 3 

3.28X10" 5 

Porous region 

11 





1007 

.61 

1.55 


2 





1004 

1.48 

3.76 



12 





1007 

.73 

1.85 


3 





2022 

1.80 

4.57 



13 





2022 

1.72 

4.37 


4 





2022 

4.25 

10.80 



14 





2022 

1.50 

3.81 


5 





4993 

3.43 

8.71 

Fractured porous region 

15 





5000 

2. 14 

5. 44 


6 





4993 

4.32 

11.00 

Fractured porous region 

16 


r 

V 

4993 

2.68 

6.81 

f 



— 

— 

- 


o 

C 






r 



! o 

c 


7 

1200 

922 

303 

1.14X10" J 

2.90xl0 -t> 

Porous region 

17 

1150 

1 894 

294i 0. 36x10 

0. 91x10 

Grain-boundary penetration 

8 





2022 

4.82 

12.20 

Fractured porous region 

18 

1150 

894 

294 

.51 

1.30 


9 





2022 

4.55 

11.60 



19 

1200 

922 

211 

.49 

1.24 


10 





4993 

10.33 

26.20 

1 


20 





290 

.56 

1.42 


11 





4993 

5.30 

13.50 



21 





305 

.63 

1.60 


t 


— 

— 

— 






22 





305 

.75 

1.91 


1 12 

1300 

977 

303 

1.36X10" 3 

3.45X10" 5 

Porous region 

23 





1005 

1.31 

3.33 


13 





303 

1.80 

4. 57 

Porous region 

24 





1007 

.92 

2.34 


14 





1004 

2.42 

6. 14 

Fractured porous region 

25 





1007 

1.09 

2.77 


15 





1004 

2.04 

5. 18 



26 





2022 

2.25 

5.72 


16 





2012 

5.82 

14. 80 



27 





2022 

2.12 

5.38 


17 





2012 

7.93 

20.10 



28 





5000 

3.01 

7.65 


18 





4622 

21.70 

55. 10 



29 


V 



4993 

3.63 

9.22 

V 

19 





4993 

18.38 

46. 68 





30 

1250 

950 

294 

0. 78X10* 3 

1. 98X10* 3 

Grain-boundary penetration 

31 

1250 

950 

294 

.83 

2.11 



32 

1300 

977 

16 

.17 

.43 



33 

1300 

977 

49 

.27 

.69 









AM 

-350 

1 

1100 

866 

295 

0. 07xl0" 3 

0. 18X10" 5 

Wall recession 

2 





295 

.27 

.69 


3 





1999 

.32* 

.81 


4 





4007 

.51 

1.30 


5 





5015 

1.90 

4.83 


6 

1150 

894 

306 

0. 39X10* 3 

0. 99X10* 5 

Wall recession 

7 

1200 

922 

295 

.29 

.74 

8 





295 

.32 

.81 


9 





301 

.41 

1.04 


10 





311 

.39 

.99 


11 





311 

.22 

.56 | 


12 





1999 

7.00 

17.80 


13 





4007 

10.40 

26.40 


14 



1 


5015 

9.80 

24. 90 ' 


15 

1250 

950 

306 

0.56X10' 3 

1.42X10* 5 

Wall recession 

16 

1300 

977 

295 

.50 

1.27 


17 





295 

.53 

1.35 



18 





311 

.90 

2.29 



19 





311 

.49 

1.24 



20 





1999 

7.60 

19.30 



21 





4007 

15.50 

39.40 



22 





5015 

16.90 

42.90 



23 





5015 

16.40 

41.70 




HS-25 


1 

1000 

811 

309 

0. 46X10" 3 

1. 17X10* 5 

Porous region 

2 

1100 

866 

294 

.85 

2. 16 



3 


294 

1.00 

2.54 



4 


294 

1.12 

2.84 



5 


990 

2.79 

7.09 



6 


1005 

3.06 

7.77 



7 


2010 

3.09 

7.85 

Fractured porous region 

8 


4002 

7.17 

18.20 

Fractured porous region 

9 


5000 

10.33 

26.24 

Fractured porous region 

10 

1200 

922 

27 

0.53X10' 3 

1. 35X10* 5 

Porous region 

11 


50 

.53 

1.35 



12 


100 

1.04 

2.64 



13 


305 

1.02 

2.59 



14 


306 

1.29 

3.28 



15 


306 

1.34 

3.40 



i6 i 


1005 

4.69 

11.90 

Fractured porous region 

17 


2010 

4.74 

12.00 



18 


2710 

9.99 

25.40 



19 


5000 

9.99 

25.40 



20 

1300 

977 

50 

0.41X10” 3 

1.04X10 -5 

Porous region 

21 





100 

.75 

1.91 

Porous region 

22 





305 

5.03 

12.80 

Fractured porous region 

23 





2010 

21.90 

55.60 



24 





4002 

33.00 

83.80 



25 





4398 

32.80 

83.30 



26 





5000 

34. 10 

86.60 





Visual analysis. - When it was removed from the capsules following a test, the 
mercury had a bright, silvery luster. Shortly after exposure to air, the mercury began 
to discolor, and a sludge appeared on the surface. The fluidity of the freshly drained 
mercury varied markedly and was dependent on material and test conditions. In all 
cases, a decrease in fluidity and an increase in sludge formation were associated with 
an increase in the severity of corrosion attack on the capsule wall in the condensing re- 
gion and an increased amount of crystalline deposit at the liquid-vapor interface. 

In general, freshly sectioned capsules showed complete "tinning” of the wall by the 
mercury. Crystalline deposits were found along the capsule walls in the liquid section, 
predominantly at the vapor - liquid pool interface. Qualitatively, there was a close 
relation between deposit buildup and corrosion penetration of the upper wall. The typi- 
cal appearance of sectioned capsules following the removal of residual mercury is shown 
in figure 3. 



Material: SICR0M0-9M 

Test temperature: 1100° F (866° K) 
Test period: 5000 hours 



AM-350 
1300° F (977° K) 
2000 hours 



H-8187 

1300° F (977° K) 
5000 hours 



067-3799 

HS-25 

1200° F (922° K) 

5000 hours 


Figure 3. - Typical appearance of sectioned capsules after removal of mercury. 


Metallographic analysis. - The maximum corrosion attack always occurred in the 
upper portion of the capsules. The variation of corrosion penetration down the capsule 
wall was wedgelike; penetration decreased in a regular manner from the top to the liquid 
pool (fig. 4). In general, examination of the corroded areas revealed no aberrations 
(i. e. , no localized region of attack with penetration much greater or much less than the 
surrounding region). The exceptions to this were the AM-350 and AM-355 alloys. In 
some instances, occasional localized deep penetrations were observed. These anom- 
alous penetrations were disregarded in estimating maximum corrosion penetration. 
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Corrosion zone 


f- 4x10 3 in. 

/ (lxlO -2 cm) 





Figure 4. - Typical variation of corrosion penetration 
down capsule wall. Photomicrographs of longitudinal 
section of capsule wall: HS-25; test temperature, 1100° F 
<866° K) ; test time, 990 hours; un etched. 
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CO 


Material: 

Test temperature: 
Test period: 



(a) Uniform wall recession. 


4xl0" 3 in. 
(Ixl0"2 cm) 



SICROMO-9M 

1100° F (866° K) 1200° F (922° K) 

1000 hours 2000 hours 

Un etched Etched 


HS-25 

1100° F (866° K) 
1000 hours 
Unetched 


(b) Channel penetration. 


(c) Porous structure. 


Figure 5 - Typical appearance of mercury-corroded areas. Photomicrographs of longitudinal section of capsule wall. (Capsule inner wall surface at 
far right of each photomicrograph.) 


The corroded areas presented one of four typical appearances: (1) uniform wall re- 
cession, approximately a 130° groove angle at surface grain boundaries (e.g. , AM-350, 
fig. 5(a)), (2) channel-like perforation, approximately a 50° groove angle at surface grain 
boundaries (e. g. , SICROMO-9M, fig. 5(b)), (3) a fine porous structure (e. g. , HS-25 and 
H-8187, fig. 5(c)), and (4) a multifractured porous structure (e. g. , HS-25, and H-8187, 
fig. 3). The groove angles given are approximate, inasmuch as grain-edge rounding and 
groove asymmetry precluded accurate angle measurements. 

All materials except AM-350 and AM-355 showed a progressive change of the form 
of the corrosion-affected zone with increasing test time and temperature. The HS-25 
and H-8187 corrosion zones progressed from a smooth wall to a porous structure and 
finally to a multifractured porous structure. The SICROMO-9M corrosion zones pro- 
gressed from a smooth wall to grain-boundary penetrations which deepened and broad- 
ened with increasing time and temperature. For the SICROMO-9M alloy, in addition to 
the gross changes in the corrosion zone, grain transformations were noted (fig. 5(b)). 
Although extensive wall recession was measured, the corrosion-affected zone of the 
AM-350 and AM-355 alloys showed little change from the initial smooth wall condition 
other than a slight surface roughening and grain transformation in a narrow region, 

O O 

1x10 inch (2. 5x10 cm) wide, adjacent to the alloy -condensate interface. 

Composition of corrosion zones, deposits, an d bulk alloys . - The compositions of 
the corrosion zones, the deposits, and the bulk alloys were determined by one or more 
of the following techniques: electron-beam microanalysis, chemical analysis, and X-ray 
diffraction. 

Typical corrosion zone compositions were determined for the alloys SICROMO-9M, 
AM-350, H-8187, and HS-25 by electron-beam microanalysis scans across the appro- 
priate region of selected samples. The SICROMO-9M analysis (fig. 6(a)) shows no effec- 
tive compositional change in the corrosion zone from that of the bulk alloy. The analysis 
of AM -350 (fig. 6(b)) shows a depletion of manganese (Mn), chromium (Cr), and nickel 
(Ni) in the region adjacent to the alloy -condensate interface. Depletion of Ni, Cr, and 
to a lesser extent cobalt (Co), with a corresponding enrichment of tungsten (W), is shown 
by the analyses of H-8187 and HS-25 (figs. 6(c) and (d)). The compositions of the corro- 
sion zones of HS-25 and H-8187, as indicated by the electron-beam microanalyses, are 
listed in table IV. Since the electron-beam analyses of these two alloys account for less 
than 85 percent of the material, other means (discussed in the following paragraphs) 
were used to establish the composition of the corrosion zones. 

X-ray diffraction patterns were obtained from the corrosion zones of selected HS-25 
and H-8187 capsules. These patterns indicate the presence of the intermetallic com- 
pounds COgW and COgW. Patterns obtained from either the as-received material or non- 
corroded areas of the test capsules do not indicate the presence of Co^W or COgW. 
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TABLE IV. - COMPOSITION OF BULK ALLOY, CORROSION ZONE 


AND CORROSION DEPOSITS 


Material 

Analysis 

Co 

Fe 

Ni 

Cr 

Mn 

W 

Mo 

N+C 






£ 




HS-25 










Bulk 

Chemical 3, 

51. 1 

2.5 

9.9 

19.8 

1.2 

15.5 

— - 

— 

Deposits 

Chemical a 

53.7 

2.7 

14.5 

27.0 

1.4 

1.7 

— - 

— 

Corrosion zone 

Electron beam* 5 

35 

1 

0 

1 


44 

— 

--- 


Mass balance 0 

44 

— 

— 

— 

— 

56 

— 

— 

H-8187 










Bulk 

Chemical 3, 

64.6 

1. 1 

1.2 

20.6 

0.6 

11.8 

0. 1 

— 

Deposits 

Chemical 3, 

69.0 

1.3 

1.4 

26.7 

.6 

2.4 

.2 

— 

Corrosion zone 

Electron beam* 5 

23 

0 

0 

2 

--- 

60 

— 

— 


Mass balance 0 

44 

— 

— 

— 


56 

— 

— 

SIC ROMO -9M 










Bulk 

Chemical 3, 

— 

90.3 

— 

8.6 

0.5 

— 

0.6 

--- 

Deposits 

Chemical 3, 

— 

87.5 

— 

7.6 

--- 


.3 

— 

AM-350 










Bulk 

Chemical a 

— 

75. 1 

4. 1 

16.6 

1. 1 

— 

2.9 

0.2 

Deposits 

Chemical 3, 

— 

69.2 

5.4 

22.5 

.8 

— 

.4 



a Analyses performed on deposits of capsules 19, 10, 28, and 23, respectively, 
for HS-25, H-8187, SICROMO-9M, and AM-350. 

^Analyses performed on corrosion zone of capsules 4, 3, 23, and 18, respec- 
tively, for HS-25, H-8187, SICROMO-9M, and AM-350. 
c Calculated from composition of bulk and deposits and assumption of conserva- 
tion of mass. 


TABLE V. - CALCULATED FLUID PARAMETERS FOR REFLUX 
CAPSULE TESTS 


Capsule bottom 

Mercury vaporiza- 

Mercury film 

Average mercury 

temperature 

tion rate 

thickness 

film velocity 

°F 

°K 

in. Vhr 

3 

00 

in. 

m 

in. /sec 

m/sec 

1000 

811 

0.34 

0. 51X10 -5 

0. 39X10" 3 

0. 99XHT 5 

0. 18 

0. 46X10' 2 

1100 

866 

.45 

.67 

.43 

1 . 1 

.22 

.56 

1200 

922 

.62 

.92 

.47 

1.2 

.27 

.69 

1300 

977 

.84 

1.4 

.54 

1.4 

.35 

.89 
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Chemical analyses were made of the bulk alloys SICROMO-9M, AM-350, AM-355, 
H-8187, and HS-25 and of deposits from selected test capsules of SICROMO-9M, AM- 
350, H-8187, and HS-25. The results are given in table IV. In all cases, except for 
SICROMO-9M, there was an enrichment of the more soluble elements, Ni, Cr, or Mn, 
in the deposits as compared with the bulk alloy. In the case of SICROMO-9M, the de- 
posits were approximately the same composition as the bulk alloy. 

A mass balance calculation was made of the weight percentage of W and Co remain- 
ing in the HS-25 and H-8187 corrosion zones by using the elemental composition data for 
the bulk alloy and the corrosion deposits. The following assumptions were made: (1) 
all material removed from the corrosion zone became deposits and (2) the amount of 
elements other than W and Co in the corrosion zone were negligible. Listed in table IV 
are the mass balance estimates of the weight percentages of W and Co in the HS-25 and 
H-8187 corrosion zones. The Co-W ratio, based on the mass balance estimates, is 2. 5 
for both materials. This ratio is consistent with the X-ray diffraction results, which 
indicated the presence of both Co 3 W and Co 2 W as the major corrosion zone constituents. 

Mercury vaporization rate and conde nsate film properties . - Table V presents the 
calculated values of the mercury vaporization rate and condensate film thickness and 
velocity for capsule temperatures of 1000°, 1100°, 1200°, and 1300° F (811°, 866°, 

922°, and 977° K). Details of the calculation are discussed in appendix A. 

With the aid of an X-ray image amplifier the course of the mercury condensation 
process was followed in a typical capsule at 1200° F (922° K) over a period of 580 hours. 
During the first 16 hours of operation dropwise condensation was noted (fig. 7(a)). 

Drops formed and grew at preferential sites along the upper portion of the capsule and 
periodically ran down the wall. After 16 hours no drops were observed; the upper por- 
tion of the capsule appeared to be of more or less uniform density as viewed through the 
X-ray image amplifier (fig. 7(b)). Mercury condensation apparently occurred uniformly 
near the top of capsule and the condensate traveled down the capsule wall in film flow. 

An estimation of the film thicknes. of the mercury condensate, based on a comparative 
optical density measurement with t known thickness of lead used as the standard, indi- 
cated the film thickness to be less . lan 1x10”^ inch (2. 5x10”^ cm). This thickness is to 
be compared with the calculated film thickness, 5x10 inch (1.27x10 cm). 


ANALYSIS 

In this investigation an attempt was made to identify, for each material tested, the 
probable rate -controlling step at any time during the progress of corrosion. The method 
of identification of rate -controlling step was based on a comparison of the experimental 
kinetic constants b r and AH r , respectively, with the theoretical kinetic constants bj 
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and AHj for each of i steps of a proposed corrosion model. Agreement of the theoret- 
ical kinetic constants of a particular step with the respective experimental kinetic con- 
stants was taken to infer that that step was rate determining. In the following presenta- 
tion, the experimental kinetic constants for each alloy were obtained by a regression 
analysis data-fitting technique. The theoretical kinetic constants associated with pro- 
posed corrosion models were derived from transport theory relations for the pertinent 
boundary conditions. 

The possibility of mechanism changeover (i. e. , change of rate -controlling step) 
during the course of the corrosion tests is also evaluated. This possibility is assessed 

(1) from the experimental results by means of data grouping and statistical testing and 

(2) from a theoretical consideration of the corrosion models. 


Experimental Kinetic Constants 


Criterion for corrosion attack. - The capsule method of corrosion testing precludes 
the use of weight change as the corrosion criterion. Therefore, in this study, maximum 
corrosion penetration was used as the sole measure of corrosion attack. Two assump- 
tions were made in applying this criterion: 

(1) Penetration is proportional to the corrosion weight change. 

(2) The test material is relatively homogeneous from the point of view of corrosion 
attack (i. e. , there are no local regions of high corrosion susceptibility such as micro- 
cracks or soluble constituent concentrations). 

The first assumption is validated, at least in the case of attack leading to uniform wall 
recession, by the principle of the conservation of mass; it also finds support in the gen- 
eral observation of increasing corrosion deposit with increasing penetration. The second 
assumption is supported by the fact that except for AM-350 and AM-355 alloys, no abnor- 
mal localized regions of penetration were found in the corroded areas of the test cap- 
sules. In the case of AM-350 and AM-355, where occasional deep localized penetrations 
were observed, the anomalous penetrations were disregarded in estimating maximum 
penetration. 

Data-fitting techniques. - Corrosion data, as a function of time and temperature, 
are customarily fitted to the empirical relations: 



(5) 


log k r = a r 



T 


( 6 ) 
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where p r is corrosion penetration (or weight loss), T is the test temperature on the 
absolute scale, t is time, and k r , a r , b r , and B r are constants. The time law (i. e. , 
the relation between p r and t at constant temperature T) is termed linear when 
b r = 1 and parabolic when b r = 1/2. The constant B r in the Arrenhius relation 
(eq. (6)) is proportional to the activation energy of the overall corrosion reaction (viz, 
AH r = 4. 58X10" 3 B r kcal/mole (19. 2 B r J/mole)). 

When the test data are fitted to equations (5) and (6) by simple linear regression 
(least-squares method), the determination of b r and B r is laden with certain difficul- 
ties. First, the precision of determination of b r and B r is generally poor because 
the dependent variable p r , as a result of uncontrolled or unknown experimental factors, 
invariably has a large random error (even though the independent variables t and T 
may have relatively small errors of measurement). Second, and perhaps as significant, 
simple linear regression analysis of the experimental data yields values for b r and B r 
which are dependent, respectively, on the values chosen for t and T (i. e. , dependent 
on the arrangement of the experimental observations, ref. ,7). 

These difficulties can be minimized or avoided to some extent by the use of multiple 
linear regression analysis (ref. 7). A single regression equation containing all the inde- 
pendent variables is used to express the dependent variable. By addition of independent 
variables to the regression equation, it is possible to account for some of the unexplained 
variation in the dependent variable obtained by simple linear regression ana/ysis. This 
resulting increase in the precision of estimate of the dependent variable p is preferred 
when analyzing data which exhibit large scatter, such as in the present case. Addition- 
ally, multiple linear regression analysis affords a more exact basis for estimating b r 
and B r than simple linear regression and also serves to summarize the results of a 
large number of observations in a single convenient statement. 

It must be emphasized that a statistical analysis of the data does not miraculously 
transform "bad" data into "good" data. The purpose in this analysis is to obtain the 
best possible estimate of the dependent variable and the kinetic constants. The uncer- 
tainty in the estimates must always be kept in mind and continually viewed with reference 
to a pertinent error statistic. 

Analysis of the test data by multiple linear regression . - The corrosion data were 
treated by the technique of multiple linear regression. The Regression Analysis Pro- 
gram with Iterative Evaluation of Relations (RAPIER) (an unpublished computer program 
written by Bert Henry, Steven M. Sidik, and R. Lambird of Lewis) was used. The esti- 
mating equation used was obtained by combining equations (5) and (6) as follows: 

B r 

log p r = a r + b r log t + — - (7) 
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TABLE VI. - RESULTS OF MULTIPLE LINEAR REGRESSION ANALYSES 


Material 

Data 

(a) 


Regression coefficients 

Standard deviation of 
regression coefficients 

Standard error 
of estimate 
of logp r , 

S log P r 

Residual error 
variance, 

s ? 

log P r 

Coefficient of 
multiple 
determination^ 

R m 

m 




Sa r 


S B 

r 

HS-25 

Total 

26 

2.5731 

0.79114 

-3951.9 

0.6142 

0.0466 

554.6 

0.162 

0.0263 

0.94 


Group 1 

15 

.5458 

. 58342 

-1712.7 

.7945 

.0896 

814.2 

.137 


.80 


Group 2 

19 

2.7311 

. 92504 

-4506.8 

.6918 

.0576 

662.0 

.145 

j, c 0200 

.94 

SIC ROMO -9M 

Total 

33 

2.4707 

0.73534 

-4185.2 

0.4332 

0.0460 

417.8 


0.0196 



Group 1 

22 

3.2780 

.98216 

-5469.8 

.5299 

.0814 

551.9 

.128 

nmm 

.89 


Group 2 

20 

.6992 

.58536 

-2159.1 

.3564 

.0338 

379.1 

.067 


.97 

H-8187 

Total 

19 

1.1873 

0.68338 

-2609.7 

0.3794 

0. 1043 

822.1 

0.199 

0.0396 

0. 76 


Group 1 

11 

-.1207 

. 54068 

-1037.4 

.2286 

.1283 

817.0 

.142 

| c . 0250 

.69 


Group 2 

14 

1.2923 

1.01128 

-3744. 5 

.3222 

.1892 

871.3 

. 169 


.77 

d AM-350 

Total 

30 

3.7425 

1.0659 

-6380.2 

0.6461 

0. 0948 

911.2 

0.247 


0.86 


a Group 1, Data exhibit initial corrosion structure; group 2, Data exhibit succeeding corrosion structure, 
b 2 

R m = Sum of squares due to multiple regression/Total sum of squares of log p r values about their mean. 

C th 

Pooled residual error variance estimate = + SSg/dfj + df 2 , where SSj is residual sum of squares for j data group and 

df. is degrees of freedom. 
d AM-350 and AM-355 data combined. 











_3 

The units arbitrarily selected for the variables were as follows: p r , inchxlO ; t, hours; 
T, °K. 

The results of the regression analyses are presented in table VI. The data were 
grouped in the following categories: total, group 1, and group 2. The latter two catego- 
ries are discussed in the next section. The 23 AM-350 and the 7 AM-355 data points 
were treated as one data set, since compositionally and structurally AM-350 and AM-355 
are almost identical. 

It is apparent from table VI that for the total data (i. e. , all the data points for a 
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given material) the values of the coefficient of multiple determination R m are large and 
range from 0. 76 to 0. 94. The large values for indicate that the functional relations 
selected for the variables were reasonably good choices, in that (r^ x 100^ percent of 
the variation in the dependent variable p can be accounted for when the data are fitted to 
equation (7). 

Taken by themselves, the values obtained from the total data for B r ranging from 

-2610 to -6380, which correspond to AH values of -12.0 to -29. 2 kilocalories per 
4 5 r 

mole (-5. 02x10 to -1. 22x10 J/mole), appear to be reasonable for this coefficient. 
However, except for AM-350, the values obtained for b r do not correspond to those 
associated with the expected processes, namely, boundary or diffusional processes 
which have linear or parabolic time-law relations, respectively. The b r values (0.79, 

0. 74, and 0. 68 for HS-25, SICROMO-9M, and H-8187) lie between the expected singular 
values of 1 or 0. 5. 

Mechanism changeover . - It is probably more than coincidental that the materials 
for which anomalous values of b r are found also exhibit marked changes in the structure 
of the corrosion zone with increasing time and temperature. Oxidation studies indicate 
that mechanism changes occur as a consequence of changes in oxide-layer thickness or 
oxide -layer adherence. Analogously, it may be expected that growth or alteration of the 
corrosion zone may change the nature of the material-transport path and, hence, result 
in a change of corrosion mechanism. If, in fact, the rate process changes during the 
course of a test, for example from linear to parabolic, then the value of b f derived 
from the total test data would assume some average value b r , where 1 > b f > 0. 5. 

The hypothesis that changes in the thickness or structure of the corrosion zone are 
accompanied by changes in corrosion mechanism was tested as follows. First, the cor- 
rosion data for each material were grouped according to the corrosion-zone criterion as 
described later in this section. Second, a multiple linear regression analysis was per- 
formed on each data group. Third, the grouped-data analyses were compared to the 
analyses previously performed on the total data. For each material, data exhibiting the 
initial corrosion structure were designated as group 1 and those exhibiting the succeed- 
ing structure as group 2. Data falling within the corrosion structure region of transition 
were included in both group 1 and group 2. The grouping criterion is as follows: 
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(1) SIC ROMO - 9M 

-3 -3 

(a) Group 1, superficial penetrations, <1x10 inch (2. 5X10 cm) 

(b) Group 2, deep penetrations, >lxl0 -3 inch (2. 5X10 3 cm) 

(2) HS-25 and H-8187 


(a) Group 1, porous structure 

(b) Group 2, multifractured porous structure 

The results of the grouped-data analyses for HS-25, SICROMO-9M, and H-8187 are given 
in table VI. 

A goodness-of-fit criterion is needed to compare the grouped-data analyses with the 
total-data analyses. The coefficient of multiple determination R^, a measure of the 
proportion of the variation in the dependent variable explained by the regression equation, 
was used to validate the choice of functional relations for the variables. This coefficient 
is a relative measure of goodness of fit. However, once the form of the estimating 
equation is found to be satisfactory, an absolute measure of goodness of fit should be 
used for comparison of analyses employing the same estimating equation. Therefore, 
the residual error variance (ref. 7), a measure of the closeness with which the estimated 
values of the dependent variable agree with the experimental values, is employed herein 
as the comparison criterion for the variously grouped data. 

2 

An inspection of the values of residual error variance s^^, shows that grouping 

the data according to corrosion structure has in all cases significantly improved the pre- 
cision in the determination of log p . Although grouping of the SICROMO-9M data into 

r -3 -3 

groups of superficial penetration (<1x10 in. (2. 5x10 cm)) and deep penetration 
(>lxl0 -3 in. (2. 5X10 -3 cm)) was arbitrary, it fortuitously proved to be a good choice. 
Trial analyses performed on the SICROMO-9M data (not shown in table VI) yielded a 
poorer fit of the data when values either slightly greater or slightly less than lxl0” 3 inch 
were used to define the division between superficial and deep penetration. 

Best estimate of regression coefficients. - The grouped data for HS-25, H-8187, 
and SIC ROMO -9M give values of b r that closely approximate the expected kinetic values 
of 1. 0 or 0. 5 and also yield reasonable values of AH r . In view of this and since data 
grouping can be rationalized on both a structural and a statistical basis, the best esti- 
mates of the coefficients of equation (7) were considered to be those from the grouped- 
data analyses. 

For AM-350 no breakdown of the data into groups was warranted. In this instance, 
for the total data the value of b r obtained was approximately 1. Further, no marked 
changes were observed in the AM-350 corrosion zone over the time and temperature 
range of the tests. The coefficients for AM-350, therefore, were considered to be those 
for the total data (table VI). 

For later comparison of results, the best-estimate values of b r and AH r and 
their standard deviations are listed in table VII. 
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TABLE Vn. - EXPERIMENTAL KINETIC CONSTANTS 


Material 

Data 

Time -law constant, 
b r 

Activation energy, 
-AH r 

kcal/mole 

J/mole 

HS-25 

Group 1 

0. 58±0. 09 

7. 8±3. 7 

33±16xl0 3 


Group 2 

. 93±0. 06 

20. 6±3.0 

86±13 

SIC ROMO -9M 

Group 1 

. 98±0. 08 

25. 0±2. 5 

105±11 


Group 2 

. 59±0. 03 

9. 9±1. 7 

41 ±7 

H-8187 

Group 1 

. 54±0. 13 

4. 8±3. 7 

20±16 


Group 2 

1.01±0. 19 

17. 1±4. 0 

72±17 

AM-350 

Total 

1. 07±0. 09 

29. 2±4. 2 

123±17 


Theoretical Kinetic Constants 

Proposed corrosion models. - Two corrosion models are presented. Each proposed 
model fits one of the two general modes of corrosion attack observed in the test capsules, 
namely, uniform attack or selective attack. In these models, movement of atoms from 
the bulk solid material to the bulk liquid mercury occurs through a series of consecutive 
and parallel material-transport steps. 

In the following presentation, S, is used to designate the i in step and S. . is used 

LL 1 1L lj K 

to designate the k parallel path for the i in step. Shown in figure 8(a) is a schematic 
representation of the transport steps proposed for the selective attack model: 

(1) Diffusion in the solid Sp consisting of the parallel steps, volume diffusion S^ ^ 

and grain-boundary diffusion g 

(2) Solution S 2 , consisting of the parallel steps, grain solution S 2 ^ and grain- 

boundary solution S 2 2 

(3) Diffusion in a static liquid Sg 

(4) Convective diffusion S^ 

Figure 8(b) illustrates the transport steps proposed for the uniform attack model (viz, 

Sfj S 2 , and S 4 ). 

General conditions and assumptions . - During a test, mercury vapor condenses con- 
tinuously near the top of the capsule (see fig. 9). The condensate flows down the capsule 
wall as a thin film in laminar steady-state flow and enters the liquid pool at y = -l . It 
was assumed that condensation at the mercury film - mercury vapor interface x = h is 
negligible. Since the measured temperature variation along the upper capsule wall is 
small and the heat flux is low, the film temperature was assumed to be uniform. 
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X 


Figure 9. - Mercury flow in reflux capsule. 



In the region of condensation the concentration of dissolved material was assumed to 
be approximately 0 (i. e. , c « 0 at y = 0). Lastly, it was assumed that steady-state 
conditions prevail at any point within the mercury film (i. e. , 3c/3t = 0). 

Time -law constant . - A mass-transport - time relation for each step of the proposed 
models was derived in the form of equation (1) from a consideration of (1) Fick’s first 
and second laws of diffusion, (2) the Navier-Stokes equation, (3) the theory of rate pro- 
cesses, and (4) the conservation of mass. The boundary conditions for each material- 
transport step are the conditions which prevail when the step under consideration is as- 
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sumed to be rate controlling. Only the overall step was considered in the derivation of 
the time -law constant bj. Parallel paths were not treated separately, since the mass- 
transfer - time relation for each path takes the same form (see appendix B). A summary 
of the estimated b. values for material-transport steps Sj, 83, S3, and S4 is given in 
table vm. 


TABLE vm. - THEORETICAL KINETIC CONSTANTS 


Material- 

Time -law 

Iron-base alloys 

Cobalt -base alloys 

transport 

step 

constant, 

b i 

Activation energy, AH i 



kcal/mole 

J/mole 

kcal/mole 

J/mole 

s i 

0.5 

30 

126 X10 3 

30 

126 X10 3 

S 2 

1.0 

11.0 

46.0 

16.4 

68.5 

S 3 

.5 

7.0 

29.3 

9.7 

37.6 

S 4 

1.0 

17.3 

72.3 

22.7 

94.9 


Diffusion in the solid S^: If diffusion through the bulk solid is the rate -controlling 
step, and with the assumption that the solid is a semi-infinite medium, the boundary 
conditions are as follows: 

for t = 0, c = c^ at x<0 and 0 > y > -l 

for t> 0, c w «0 at x = 0 and 0 > y > -l 

where c^ and c w are the concentrations in the bulk and at the wall, respectively, of 
the diffusing species. The solution of this problem is 

= “l c b D s /2 ' 1/2 <*> 

where is a constant, D g is the diffusion coefficient for the species in the bulk solid 
(assumed to be independent of c). For the selective attack model, if it is assumed that 
the penetration is proportional to w^ equation (8) may be written as 

Pi = ^D^t 1 / 2 (9) 

where Pi is the penetration and ( 3 ^ is the appropriate proportionality constant. For 
the uniform attack model, if it is assumed that recession of the wall is proportional to 
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the amount of diffusing substance, solution of the problem of a moving boundary (refs. 8 
and 9) leads to a relation identical to equation (9) with the exception that the proportion- 
ality constant is different and p^ represents the amount of wall recession. 

Solution Sg: E solution is the rate -controlling step, the boundary conditions for 
t > 0 are c = c(y) at x = 0 and 0 > y > -l. The solute concentration c(y) is a function 
of y as determined by the convective dEfusion parameters of the flowing mercury film. 
The amount of a given species that passes into solution per unit time per unit surface 
area of solid is 


G = k sr [ c o " c( y>] (10 > 

where k is a solution rate constant (a function only of temperature) and c is the 
sr u 

equilibrium solubility concentration of the solute in mercury. 

H the rate of corrosion penetration or wall recession is assumed to be proportional 
to the flux of atoms through the surface into the mercury, then the surface flux is given 
by 


dp 2 

3x = o = “ 2 ir 


(ii) 


where is a proportionality constant and p 2 is the corrosion penetration or wall re- 
cession. For steady-state conditions, 

j x=0 = G 

Equating equations (10) and (11) and integrating with respect to t yield the following 
kinetic relation for the solution-controlled process: 


P2 = — [ c o " 

a 2 


( 12 ) 


DEfusion in static liquid S^: E dEfusion along quiescent liquid-filled channels or 
networks is the rate -controlling step, the boundary conditions for t > 0 are: 

c = c Q at x = pg and 0 > y > -l 

c = c’(y) at x = 0 and 0 > y > -l 
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where p 3 is the depth of greatest penetration of the channel and c'(y) is the solute con- 
centration at the channel mouth (a function of y as determined by the convective diffu- 
sion parameters of the mercury film). 

The flux of solute from the head to the mouth of the channel is inversely proportional 
to the length of the diffusion path and directly proportional to the concentration difference 
as follows: 


D f[ c o - c, (y)] 

P3 


(13) 


where is the diffusion coefficient for the solute in mercury assumed to be indepen- 
dent of c . 

If the penetration rate is assumed to be proportional to the flux, then 


j = a 3 


dpg 

dt 


(14) 


where oig is a proportionality constant. Substituting the relation for j of equation (13) 
into equation (14) and integrating with respect to t result in the kinetic relation for the 
static liquid-diffusion-controlled process 

|2D,[c 0 - c'(y)] 

P3 t ^ 

Convective diffusion S^: If convective diffusion is the rate -controlling step, the 
boundary conditions for t > 0 are, effectively, 

c = c q at x = 0 and 0 > y > -l 

and 



— = 0 at x = h and 0 > y > -l 
dx 

The concentration distribution for a solute diffusing through a flowing thin film in 
laminar flow and steady state is given in reference 10 as 
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( 16 ) 


3vf—~ — \ i£ = D, — 

\ h 2h 2 / 9 y ax 2 

where v is the average film velocity (eq. (A3)) and h is the film thickness (eq. (A4)). 

Since the liquid film is flowing at relatively low velocities (see table V), a fully developed 

diffusion regime will be established across the thickness of the film a short distance from 

the top of the film. Calculations indicate that the thickness of the diffusion boundary layer 

-3 -2 

becomes comparable to the film thickness at y « 4x10 inch (10 cm). 

The convective diffusion mass transfer problem of S 4 is analogous to the problem 
of heat transfer to a fluid in laminar flow along a wall which is maintained at uniform 
temperature. The solution to the heat-transfer problem (ref. 11), when applied to the 
convective diffusion problem, yields the following expression for the flux of mass at the 
surface : 


j 


x=0 


D 7 Cf 

l o 
h 


e" n ^ 


(17) 


_n 

for values of | greater than 0. 12 (i. e. , y > 10 cm) where n and f are constants 
and £ = D^y/3vh 2 . 

If it is assumed that the penetration rate is proportional to the flux, 


■*x=0 


= «/ 


dp 4 

dt 


(18) 


where a 4 is a proportionality constant and p 4 is the penetration. Substituting the ex- 
pression for the flux (eq. (17)) into equation (18) and integrating with respect to time 
yield the kinetic relation for the convective diffusion controlled process 
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(19) 


Activation energy . - The activation energy AH^ associated with each material- 
transport step is defined by equation (2). For purposes of evaluating AHj, the differen- 
tial form of equation (2), namely, 


R d(ln kj) 


= AH. 


d — 


(4) 
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was used. Expressions for k^, kg, kg, and k^ (multiplied by a proportionality constant) 
are obtained directly from equations (9), (12), (15), and (19), respectively. The expres- 
sions for kp evaluated at y « 0, were introduced into equation (4) and the differential 
operation was performed. The AH^ values obtained are as follows (see appendix C for 
details of the calculation): 


-AH X = 0.5Qj 

(20) 

iH 2 - 4 + 4r 

(21) 

= 0. 5 (q* + Q y + RT) 

(22) 

- Qs + Q v + 2 - 75 RT 

(23) 


where Q , is the effective activation energy for chemical diffusivity in alloy A (where 

a A A 

A represents any one of the four alloys tested), Q„ and Q ct , are, respectively, the 

effective activation energy of solution and solution rate in mercury of alloy A, and 

is the activation energy for viscous flow in mercury. The average temperature over the 

experimental range, 900° K, is used for T. 

Of the four activation energies only can be readily estimated. The value of 

was assumed to be 1. 3 kilocalories per mole (5. 4x10^ J/mole) (ref. 12). Because of the 

complex nature of the diffusion and solution processes of multicomponent systems, sev- 

A A A 

eral simplifying assumptions must be made to obtain estimates of Q^, Q sr > and Q g . 

As a result, too much confidence should not be placed in the value calculated for AH,. 

A 1 

The evaluation of was based on the following assumptions: 

(1) Volume diffusion ^ is the predominant mode of diffusion under the test con- 
ditions. 

(2) The test alloys can be treated as pseudobinary systems (i. e. , Co-Cr or Fe-Cr 
systems). 

(3) The thermodynamic activity coefficients of the major constituents of each alloy 
are effectively constant with change in composition. 

The activation energy for volume self -diffusion of the major constituent elements of the 
test alloys (Fe, Co, Cr) have comparable values (ref. 13) that fall within the range of 
60±10 kilocalories per mole (2. 51±0. 42x10^ J/mole). Therefore, the value of the effec- 
tive activation energy of chemical diffusivity for all the test alloys was taken to be ap- 
proximately 60 kilocalories per mole (2. 51x10^ J/mole). 

A 

The evaluation of Q_ was based on the following assumptions: 
s 
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(1) Only the base element of any given alloy (i. e. , Fe or Co) need be considered in 
the solution process. 

(2) The influence of mutual solubility effects on the activation energy of solution is 
negligible. 

The first assumption is supported by the fact that the alloy base element was the major 
component of the corrosion deposits (see table IV). The second assumption is supported 
by the results of a study of the solubility of the alloys HS-25, H-8187, and SICROMO-9M 
in mercury (ref. 14). 

The data from reference 14 for the cobalt solubility of the alloys HS-25 and H-8187 
as a function of temperature are presented in figure 10. Plotted on this same figure are 
the data from references 15 and 16 for elemental cobalt solubility. The data from refer- 
ence 14 for iron solubility of the alloy SICROMO-9M are shown in figure 11. Also plot- 
ted on figure 11 are the data from references 15 and 17 for elemental iron solubility. 
Although AM -350 data are given in reference 14, these data are not presented herein be- 



Inverse temperature, “K"- 1 

Figure 10. - Solubility of cobalt in mercury. For elemental cobalt, log (ppm Co) = 3. 8725 
- (3585. 3/T); standard deviation of log (ppm Co), 0. 1925; standard deviation of slope, 560.8. 
For cobalt from HS-25, log (ppm Co) = 4. 1730 - (4093. 4/D; standard deviation of log {ppm Co), 
0.2029; standard deviation of slope, 1099. 
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Figure 11. - Solubility of iron in mercury. For elemental iron, log (ppm Fe) 

= 2.4396 - (2409. 1/T); standard deviation of log (ppm Fe), 0. 1911; standard devi- 
ation of slope, 453.9. For iron from SICROMO-9M, log (ppm Fe) = 0.9442 
- (1426. 1/T); standard deviation of log (ppm Fe), 0.2826; standard deviation of 
slope, 752. 


cause of their extremely large scatter. Our estimates of the regression line and rele- 
vant standard deviations for the alloy and element data are shown on each figure. Since 
there are only four data points for the H-8187 alloy, no regression estimate was made. 
The slopes of the regression lines shown in figures 10 and 11, when multiplied by the 
factor 4. 58, yield the appropriate activation energy of solution in units of kilocalories 
per mole. 

The hypothesis that the slope of the solubility line for the base element of an alloy is 
the same as the slope of the solubility line for the pure element (i. e. , the second as- 
sumption) can be tested by analysis of covariance (ref. 18). The results of a comparison 
of the slopes of the regression lines are as follows. For the HS-25 alloy slope compared 
to the elemental cobalt slope, the 95-percent confidence interval for the difference be- 
tween slopes is 0. 51±2. 53xl0 3 . For the SICROMO-9M alloy slope compared to the ele- 
mental iron slope, the 95 -percent confidence interval for the difference between slopes is 
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3 

0. 98±1. 86x10 . Since the interval in both cases includes zero, the hypothesis cannot be 
rejected. 

For the reasons just discussed, the activation energy of solution for the cobalt- and 
iron -base alloys was taken to be 

qHS-25 w qH-8187 w q Co _ 16 4 kcal / mole (68. 5x io 3 j/mole) 
s s s 

q! M » qAM- 350 « Qj e = 11.0 kcal/mole (46. OxlO 3 J/mole) 
s s s 


No data are available for solution rate activation energy in mercury. This activa- 
tion energy is probably related to the energy needed by the dissolving species to diffuse 
through an impurity or reaction-product barrier film at the solid-liquid interface. In 
the tests reported herein, there was no deliberate addition of a mercury solution rate 
inhibitor (e. g. , zirconium or titanium, which tend to form nitride surface films (ref. 19)). 
Further, X-ray observations showed good wetting of surfaces in contact with mercury, 

which suggests the absence of an oxide or oil impurity film at the solid-mercury inter- 

A 

face. It was therefore assumed that Q is zero. 

far 

Listed in table VIII are the AH^ values for material-transport steps S^, S 2 , Sg, 
and calculated from equations (20) to (23) and the appropriate values. 

Variation of penetration as function of y. - Where the concentration of dissolved 
material is also a function of y (i. e. , rate control by S 2 , Sg, or S^), a variation of 
corrosion penetration down the capsule wall may be expected. Corrosion penetration 
should decrease in a regular manner (determined by the solution rate constants and con- 
vective diffusion parameters) from y = 0 where c(y) = 0 to y = -l where c(y) < c . 

Mechanism changes. - The proposed corrosion model for selective attack implicitly 
presupposes a change of mechanism during the progress of corrosion. Phenomenologi- 
cally, this can be viewed as a consequence of changes wrought in the length and nature of 
the material-transport path. Analytically, this is made evident when the kinetic relation 
for the overall corrosion reaction (i. e. , for all steps of the model viewed together) is 
examined. The overall kinetic relation for the selective attack model is 

at = y- jP 2 + y 2 P + r 3 P 2 + r 4 p (24) 


where y. , y^, y g, and y ^ are constants related to the respective constants of the 
transport steps, as described in the section Time-law constant, and p is the corrosion 
penetration predicted by the selective attack model. (See appendix B for derivative of 
eq. (24). ) It is seen from equation (24) that when p is small either S 2 or S4 (depend- 
ing on the relative values of yg an( * 7 4) must be the rate -determining step. Depending 
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on the relative values of the constants, either S-^ or Sg will become rate controlling as 
corrosion penetration increases with time. To recapitulate, the selective attack model 
predicts that the initial rate -controlling step follows a linear corrosion rate (associated 
with a boundary-region process). Further, during the progress of corrosion, this initial 
rate-controlling step is succeeded eventually by a step that follows a parabolic corrosion 
rate (associated with a diffusion process). 

No material-transport path change is postulated for the uniform attack model. 
Therefore, where the uniform attack model applies, no mechanism change may be ex- 
pected. 


DISCUSSION OF RESULTS 
Corrosion Rates and Rate-Controlling Steps 

Inferences of the rate and the rate -controlling step in the corrosion attack of mer- 
cury on the materials tested are presented herein. These inferences are based primar- 
ily on a comparison of the experimental and theoretical kinetic constants that are derived 
in the ANALYSIS section and listed in tables VII and VIII, respectively. Subsidiary evi- 
dence as to the nature of the rate -controlling step is found in a comparison of the ob- 
served and predicted variation of corrosion penetration along the capsule wall. The uni- 
versally observed wedge -like penetration pattern in longitudinal sections of corroded cap- 
sules corresponds to the predicted variation of corrosion penetration when any step but 
is rate controlling. 

Cobalt -base alloys, HS-25 and H-8187 . - In the initially observed corrosion regime, 
data group 1 , the corrosion rate was parabolic, and the probable rate -controlling step 
was Sg. In the succeeding corrosion regime, data group 2, the corrosion rate was lin- 
ear. In the case of both alloys, the AH values fell between the AH. values for S 0 

r 1 z 

and S^. Because of the relatively small difference between AH. for S 2 and S 4 , the 
uncertainty in the estimate of AH r , and the assumption made in estimating the AH^ 
values, an unequivocal choice of rate -controlling step cannot be made. 

Iron-base alloy, SICROMO-9M . - In the initially observed corrosion regime, data 
group 1, the corrosion rate was linear. The probable rate -controlling step was taken to 
be S^. The judgment was somewhat arbitrary because, although the AH r value was 
closer to AH^ than to AHg, it was significantly larger than AH^. (See the discussion 
for AM-350. ) In the succeeding corrosion regime, data group 2, the corrosion rate was 
parabolic, and the probable rate -controlling step was Sg. 

Iron-base alloy, AM-350 . - Over the entire test range, the corrosion rate was lin- 
ear. The probable rate -controlling step was S^. Again, as in the case of SICROMO-9M, 
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this judgment is somewhat arbitrary. The result for AM-350 and SICROMO-9M 

(AH r > AH^ by 8 to 12 kcal/mole (3xl0 4 to 5x10^ J/mole)) is similar to that reported by 

Horsley (ref. 4) for the corrosion of mild steel in mercury from 900° to 1200° F (755° 

to 922° K), assuming boundary-layer control. From Horsley's calculations AH was 

4 ^ 

greater than AH^ by 11 kilocalories per mole (4. 6x10 J/mole). These discrepancies 
between experimental and theoretical activation energy may possibly reflect inaccuracy 
in the iron solubility data or, conversely, an incorrect corrosion model. 


Mechanism Change and Corrosion Regimes 

Analysis of the experimental data indicates that mechanism changes occur during 
the progress of corrosion for all alloys studied, with the exception of AM-350. Theoret- 
ical consideration of the corrosion models support these results. 

When a change in mechanism takes place, the corrosion reaction enters a new cor- 
rosion regime. Three corrosion regimes can be characterized from this present study: 
linear I, parabolic, and linear II. In the following exposition of the corrosion regimes, 
figure 12 is used as a graphic example. In this figure are the 1100° F (866° K) corrosion 
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rate curves, derived from equation (7), and the pertinent values found in table VI for 
HS-25, H-8187, SICROMO-9M, and AM-350. The solid portion of each curve covers the 
test range. The dashed portions are conjectural. 

Linear I . - At the start of each test and for some time thereafter the inner surface 
of a capsule is relatively smooth, with no deep corrosion penetrations. As such, the 
corrosion reaction is best represented by the uniform attack model. Accordingly, the 
initial corrosion rate is expected to be linear for all materials. 

A linear rate is maintained when the capsule wall exhibits a continuous uniform re- 
cession. Such behavior is exemplified by the AM-350 alloy, as indicated in figure 5. 
Attack on the AM-350 alloy progresses by grain -boundary penetration. Since, however, 
the groove angle formed is approximately 130°, the corrosion attack moving on a wide- 
angle front results effectively in wall recession. 

Parabolic . - Continuing corrosion attack on the materials tested, other than AM-350, 
results in the formation of a channel-penetrated or porous-network type of corrosion 
zone. This developing corrosion zone interposes itself as a new link in the mass- 
transport chain. Transport of matter across the channeled or porous region occurs by 
diffusion along the quiescent liquid paths created by mercury contained within the chan- 
nels or the network. When the corrosion channels or network reach a sufficient depth, 
the rate of material transport along the liquid diffusional path becomes slower than along 
other paths in the transport chain. Thus, the rate -controlling step changes, and the ini- 
tially observed linear corrosion rate becomes parabolic. The change of rate from linear 
to parabolic is continuous since the physical changes occurring in the corrosion zone are 
continuous and progressive. For simplicity of presentation, no attempt was made to in- 
dicate a continuous change in rate on the pertinent curves in figure 12. 

The formation of channels in the SICROMO-9M alloy can be attributed to intergranu- 
lar attack advancing with a small groove angle, approximately 50°. No satisfactory 
explanation can be offered for the pronounced difference in magnitude of the groove 
angles produced in the SICROMO-9M and AM-350 alloys, which in the case of the latter 
alloy leads to wall recession. 

The formation of a porous Co^W and COgW network in the HS-25 and H-8187 alloys 
can be attributed to leaching of the more soluble elements from the corrosion zone. 

Since Co^W and CogW are not found in the as-received material nor in the material out- 
side the corrosion zone of the tested capsules, it can be concluded that these compounds 
are formed in situ by reaction of the residual elements following the leaching process. 

Figure 12 indicates that at 1100° F (866° K) SICROMO-9M exhibits a linear corro- 
sion rate until about 1400 hours, after which the rate becomes parabolic. A comparable 
linear -to-parabolic rate change for HS-25 and H-8187 remains unconfirmed because of 
the lack of sufficient short-time data. The initial linear rate curve shown for these al- 
loys in figure 12 is conjectural. 
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Linear n. - In the case of the HS-25 and H-8187 alloys, a change in density probably 
occurs within the corrosion zone as a result of leaching. From mass balance considera- 
tions (see table IV), it can be estimated that the density of the fully developed HS-25 and 
H-8187 corrosion zones is, respectively, 0. 24 and 0. 18 the density of the original mate- 
rial. The density changes and the accompanying volume changes probably induce 
stresses within the corrosion-zone network and between the network and the base mate- 
rial. When the stresses exceed the critical shear stress, cracks develop (see fig. 3, 
p. 10). The cracks allow the flowing mercury direct access to the base material. Asa 
result, the diffusional path through the porous corrosion zone is bypassed, and the corro- 
sion rate reverts from parabolic to linear. The cracks, once initiated, can be expected 
to grow as corrosion continues and to be accompanied by spalling; thus, the linear corro- 
sion rate is maintained. This particular corrosion behavior is closely analogous to cer- 
tain oxidation reactions where parabolic-to-linear rate changes occur when developed 
stresses cause cracking and scaling of oxide layers. 

In figure 12, the corrosion rate of HS-25 and H-8187 at 1100° F (866° K) changes 
from parabolic to linear at 1100 and 2400 hours, respectively. This change is consistent 
with the first appearance of cracks in the corrosion zone of test capsules. Cracks ap- 
pear between 1000 and 2000 hours for HS-25 and between 2000 and 5000 hours for H-8187. 


Liquid Flow Velocity 

Liquid flow velocity cannot be treated as an independent variable in reflux capsule 
corrosion tests. In these tests the mercury film velocity is a function of the vaporization 
rate, which in turn is a function of temperature. Since, however, the temperature- 
velocity functional relation is relatively weak and the temperature range of investigation 
is relatively narrow, the effect of velocity on the test results was assumed to be negli- 
gible. This was justified as follows: The data from the linear rate group for each mate- 
rial were fitted by multiple regression analysis to the estimating equation 

B 

log P r = a r + b r !og t + — L - 0. 75 log T (25) 

The estimating equation was derived from equation (19) with the appropriate substitution 
for the constants D^, c Q , and h in terms of T (see appendix C). Equation (25) reflects 
the effect of velocity on corrosion with the assumption of control. The residual 
error variance and the coefficients b r and B r (eq. (25)) did not differ significantly 
from the values found in table VI for the pertinent data fit to equation (7). 

Although velocity effects were not discernible in this investigation, such effects 
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should be manifest at high liquid flow velocity. Speculations and observations on the ef- 
fect of high flow velocity on erosion and corrosion are presented in the following para- 
graphs. 

Erosion . - Flow velocities within the reflux capsules were 0. 2 to 0. 4 inch per sec- 
ond (0. 5 to 1. 0 cm/sec), too low to incur erosional effects. However, at higher velocity, 
materials inclined to intergranular penetration are likely to suffer erosion through grain 
removal. Further, as a consequence of this erosion, the mode of corrosion attack will 
differ from that observed in the capsules. In the case of SICROMO-9M, for example, 
initial intergranular attack would produce a layer of peripherally corroded grains at the 
alloy- liquid interface. Through the action of a high stream velocity these grains would 
be readily swept away. As a consequence of continuing erosional grain removal, channel 
formation would be curtailed. With channel formation curtailed, mechanism change is 
precluded, and corrosion attack would continue in the linear I regime. 

This speculation is supported by the results of SICROMO-9M forced-convection- 
boiling mercury loop tests run at liquid flow velocities as great as 108 inches per second 
(274 cm/sec) (refs. 20 and 21). No channeled-corrosion zone was found, although inter- 
granular attack was evident. Also, the wall thickness loss was considerably greater than 
that predicted by the capsule tests for comparable time and temperature conditions. 

Corrosion . - The effect of liquid flow velocity on corrosion is related to the nature 
of the rate -controlling step of the corrosion reaction. When boundary- layer diffusion 
S 4 is rate controlling, corrosion will be proportional to v 11 , where n is 1/3 for laminar 
flow and n is 1/2 to 1 for turbulent flow. When solution Sg is rate controlling, equa- 
tion (12) indicates that penetration is proportional to c Q - c(y). The term c(y) is depen- 
dent on solution rate and stream velocity (as defined by the convective diffusion param- 
eters of the stream) and can range in value from 0 to c Q . Consequently, Pg will vary 
from a maximum when c(y) = 0 to 0 when c(y) = c Q . For two cases of interest, Pg is 
independent of v: for reflux capsules at the point of maximum penetration (i. e. , y = 0, 
c(y) = 0) and for high flow, low solubility rate situations (i. e. , c(y) = constant). 

The results of this investigation suggest that, for AM-350 and SICROMO-9M alloys, 
S 4 control is likely in the linear corrosion regimes; no resolution could be reached on 
the question of Sg or S 4 control for the HS-25 and H-8187 alloys. Evidence supporting 
S 4 control is found in the results of tests of forced-convection corrosion loops fabricated 
of SICROMO-9M alloy (refs. 20 and 21) and HS-25 alloy (unpublished results of tests by 
A. Vary of Lewis). In all these tests, mercury saturation temperature was about 
1100° F (866° K). The velocity range covered was from 12 to 108 inches per second 
(30. 5 to 274 cm/sec). Although the penetration data as a function of velocity exhibit a 
large variability, in all cases corrosion penetration is greater than that expected for 
linear corrosion rates with Sg control. 
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SUMMARY OF RESULTS 


A total of 108 mercury reflux capsule tests were performed on the cobalt -base al- 
loys HS-25 and H-8187 and the iron-base alloys SICROMO-9M, AM-350, and AM-355. 
(The latter two alloys are treated as one alloy. ) The test temperature range was 1000° 
to 1300° F (811° to 977° K), and test times were as long as 5000 hours. 

Metallographic analysis revealed that corrosion resulted in either uniform attack, 
characterized by wall recession (AM-350), or selective attack, characterized by a 
leached-out porous corrosion zone (HS-25 and H-8187) or a channel-like grain-boundary- 
penetrated corrosion zone (SICROMO-9M). 

Values of the experimental kinetic constants (i. e. , the time-law constant b p and the 
activation energy AH r ) were obtained from a multiple regression fit of the corrosion 
data to the estimating equation 


AH 

log P = a + b log t + 

r r r 2. 303 RT 

where p f is the observed corrosion penetration, t is test time, T is test temperature, 
R is the gas constant, and a r is a constant. 

The theoretical kinetic constants b^ and AH^, for each of i steps of proposed 
corrosion models, were obtained from transport theory relations and the pertinent 
boundary conditions. 

Analysis of the experimental data indicated that mechanism changes occurred during 
the progress of corrosion for all alloys studied, with the exception of AM-350. Theo- 
retical consideration of the corrosion models supported these results. 

Inferences of rate and rate -controlling step were made by comparing AH^ and b^ 
for each step in the proposed corrosion model with AH r and b r , respectively. Because 
of the lack of data needed to estimate the activation energy for several pertinent pro- 
cesses, there is some uncertainty in the estimated AHj^ values. Consequently, the in- 
ference of rate -controlling step must be considered as tentative. 

Two corrosion regimes were inferred from the HS-25 and H-8187 alloy data. The 
first was a parabolic rate regime with liquid diffusion as the probable rate -controlling 
step. This regime was succeeded by a linear rate regime with boundary -layer diffusion 
as the probable rate -controlling step. Although not experimentally verified, there is 
good reason to believe that the parabolic regime is preceded by a linear rate regime. 

Two corrosion regimes were inferred from the SICROMO-9M alloy data. The first 
was a linear rate regime with boundary-layer diffusion as the probable rate-controlling 


38 



step. This regime was succeeded by a parabolic rate regime with liquid diffusion as the 
probable rate-determining step. 

A single corrosion regime of linear rate with boundary -layer diffusion as the prob- 
able rate -controlling step was inferred in the case of the AM-350 alloy. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 24, 1967, 

120-27-02-18-22. 
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APPENDIX A 


MERCURY VAPORIZATION RATE AND CONDENSATE FILM PROPERTIES 

A determination was made of the temperature profile along the capsule and the heat 
loss from the top of the capsule at several temperature levels in order to estimate the 
mercury vaporization rate in the reflux capsule tests. A typical capsule was used with 
five thermocouples spot welded to the wall in the positions shown in figure 13. Other 
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Distance from bottom of furnace, cm 

Figure 13. - Typical temperature profiles for a capsule containing mercury. 
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Thermocouple position 



1 2 3 4 5 6 

Distance from bottom of furnace, cm 

Figure 14. - Typical temperature profiles obtained with an empty 
capsule. 
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thermocouples were distributed throughout the capsule furnace assembly. The setup was 
operated at 1000°, 1100°, 1200°, and 1300° F (811°, 866°, 922°, and 977° K) with the 
capsule one-third filled with mercury. The capsule was then removed, emptied of mer- 
cury, and resealed in vacuum with an electron-beam welder. The capsule was re- 
installed in the furnace and operated over a range of temperature from 900° to 1600° F 
(755° to 1144° K). 

The temperature profiles for the runs made with mercury in the capsule are shown 
in figure 13. The slight temperature rise indicated along the lower two -thirds of the 
capsule is probably a result of the proximity of thermocouples 2 and 3 to the heater. The 
temperature drop observed over the top one-third of the capsule is most likely somewhat 
larger than the true drop, because of heat-conduction loss from convective cooling of the 
thermocouple wire. 

Temperature profiles obtained with the evacuated empty capsule are shown in fig- 
ure 14. The profiles show a nearly uniform temperature for the lower portion of the 
capsule and a temperature drop of 200° to 400° F (111 0 to 222° K) for the upper portion 
of the capsule. The anomalous temperature drop indicated by the top thermocouple (5) 
probably results from two factors. The first factor is an added resistance to heat con- 
duction at the thermocouple location produced by the dimple machined in the surface of 
the capsule. The wall thickness at the dimple was 0. 010 to 0. 020 inch (0. 025 to 
0. 051 cm) as compared with 0. 040 inch (0. 10 cm) for the nominal wall thickness. The 
second factor is an enhanced heat-conduction loss down the thermocouple wires caused 
by convective cooling. It is assumed that the true top temperature is that value obtained 
by straight-line extrapolation of the temperatures given by thermocouples 3 and 4. 

The heat loss from the top of the capsule was calculated from the following equation: 

q = kA £1 (Al) 

L 

where k is the thermal conductivity, A is the cross-sectional area of the capsule wall, 
and AT/L is the slope of the linear portion of the temperature -profile curve between 
the thermocouples 3 and 5 of figure 14. Heat loss from the capsule top is plotted in fig- 
ure 15 as a function of true top- surface temperature of the empty capsule. 

Since the heat loss from the top of the capsule for any given top temperature is inde- 
pendent of whether the capsule is filled or empty, values of heat loss for the filled cap- 
sule tests can be obtained directly from figure 15. 
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Figure 15. - Heat loss from capsule top as function of temperature. 


The vaporization rate m was obtained from the relation 

m=-SL. 

V 


(A2) 


where H v is the latent heat of vaporization and p is the density of mercury. Calculated 
values of m are given in table V for capsule bottom temperatures of 1000°, 1100°, 

1200°, and 1300° F (811°, 866°, 922°, and 977° K). 

Average velocity v over the film cross section is (ref. 11, p. 671) 



(A3) 


where h is the film thickness, v is the kinematic viscosity, and g is the gravitational 
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constant. A relation for the film thickness can be obtained as follows. The average 
velocity is equal to the vaporization rate m divided by the flow cross-sectional area of 
the film. Since h « d, where d is the capsule inside diameter, the flow cross- 
sectional area is approximately equal to rrdh. Substituting m/Vdh for v in equation 
(A3) and solving for h yield 


h = (IS!?* 

\TTdg I 


(A4) 


The calculated film thickness and average velocity are given in table V for capsule 
temperatures of 1000°, 1100°, 1200°, and 1300° F (811°, 866°, 922°, and 977° K). 
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APPENDIX B 


DERIVATION OF KINETIC RELATION FOR SELECTIVE ATTACK MODEL 

The kinetic relation for the selective attack model can be obtained by considering the 
overall corrosion reaction in terms of a flux of matter j which has to overcome a series 
of resistances r j. Each resistance is associated with a corresponding step S^. The 
sum total of all the steps compose the mass-transfer path. The driving force for each 
step is a concentration difference Ac^. Thus, analogous to current flow in an electric 
series circuit, 


j = 


E Ac i 

E r i 


(Bl) 


When parallel paths exist, as for example in the case of volume and grain-boundary 
diffusional paths ^ and ^ i n the s °lid diffusion step (fig. 8), the total resis- 
tance r ^ is given by the expression 


r l 


1 


+ 


1 


r l,l r l,2 


Consequently, when r^ ^ ~ 2’ r l ~ r l 2^ or w ^ en r i l ^ r i 2’ r l ~ r l 2' Thus, 

the total resistance may be expressed in terms of the smaller of the parallel path resis- 
tances. 

Two types of resistances may be considered: (1) resistance r^ ^ associated with 
steps involving a boundary region (e. g. , S 2 and S^) and (2) resistance r^ ^ associated 
with steps involving diffusional processes (e.g. , and Sg). The relation between r. ^ 
and the depth of corrosion penetration p is 


r 


i , 1 


r t P 


(B2) 


i.U 

where y . is a characteristic constant for the i Ln step, 
pressed simply as 


Resistance 


rj ^ may be ex- 


r i,b = ^i 


(B3) 


When the flux of matter is constant, 
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jxpSE 

dt 


where cp is a proportionality constant. Combining equations (Bl) and (B4) and making 
the appropriate substitution for the r i terms yield the following corrosion-rate expres- 
sion for the selective attack model: 


E ac i 

= hi 

dt yjP + r 2 + r 3 P + r 4 


(B5) 


If it is assumed that Ac^ is constant, integrating equation (B5) and combining constants 
into a new constant a (the integration constant being 0 since p = 0 at t = 0) yield the 
following kinetic relation: 

2 2 

at = yjP + r 2 P + r 3 p + r 4 P (B6) 
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APPENDIX C 


ACTIVATION -ENERGY CALCULATIONS 

The following equation was used to calculate the activation energies for each 
material-transport step: 


AH i = 


R d(ln kj) 


dl 


(4) 


Listed below are the pertinent expressions for kj obtained from equations (9), (12), (15), 
and (19) and evaluated at the point of maximum attack. An appropriate proportionality 
constant is included: 


k = e 
K i e rs 

(Cl) 

k 2 = e 2 k sr c o 

(C2) 

k 3 = e 3^ D i c o^ 1/2 

(C3) 

k -e D ^° 
44 h 

(C4) 

Diffusion in Solid S^ 


1 dl 

T 

(C5) 


The temperature dependence of D is as follows: 

s 



where 
loy A. 
yield 


A 

D q is a constant and is the effective activation energy for diffusivity in al- 

Substituting for D in equation (C5) and performing the differential operation 
s 


-AHj = 0. 5 


(C7) 


Solution S2 

ah - R d ( ln £ 2 k sr c o) 

2 di 

T 


(C8) 


The temperature dependence of k and c is given by the following relations: 

sr o 


k 


sr 



(C9) 


c 


o 



(CIO) 


A A 

where X and X are constants and Q and Q are, respectively, the effective acti- 

b bl 

vation energy of the solution and the solution rate in mercury of the elements of alloy A. 

Substituting for k and c in equation (C8) and performing the differential operation 
sr o 

give 


-ah 2 = qA + q a 


(oil) 


Diffusion in Static Liquid S3 

_ R d in [g 3 (D^ Cp) 1/2 ] 


AH 3 = 


di 


(C12) 


The temperature dependence of can be approximated by the Stokes-Einstein relation: 
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(C13) 


D 


l 


V 

6l Vo 



where is the Boltzmann constant, rj Q is a constant, r Q is the radius of the diffusing 
substance (approximately the same for all metal atoms considered in this study), and 
is the activation energy for viscous flow. Combining equations (C12), (CIO), and (C13) 
and differentiating yield 


-AH 3 = 0.5(q^ + Q v + Rt) 


(Cl 4) 


Convective Diffusion S4 

Rd ( ln£4 nr) 

AH, = A ^—L (C 15) 

di 

T 

An analysis of the mercury vaporization data (table V) indicates that h is proportional 
to T*" . Substituting the appropriate temperature -dependent relations for D^, c q , 

and h into equation (C15) and differentiating give 

-AH 4 = + (^ + 2.75 RT (C16) 
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APPENDIX D 


SYMBOLS 


A cross-sectional area of capsule 
wall, 

a constant 

B r regression coefficient 

b time -law constant 

Q 

c concentration, kg/m 

c equilibrium solubility concentra- 

0 3 

tion, kg/m 

D, diffusion coefficient in liquid, 

L n 

ni /hr 

D diffusion coefficient in solid, 

s 2 

m /hr 

d capsule inside diameter, m 

f constant 

2 

G solution mass flux, kg/(m )(hr) 

2 2 

g gravitational constant, (N-m )/kg 

H y latent heat of vaporization, J/kg 

AH activation energy, J/mole 

h condensate film thickness, m 

j mass flux, kg/(m )(hr) 

k specific rate constant 

k^ Boltzmann constant, J/°K 

k solution rate constant, m/hr 

sr 

L linear portion of temperature pro- 

file, m 

l vertical length of condensate 

film, m 

3 

m vaporization rate, m /hr 


n constant 

p corrosion penetration or wall re- 

cession, m 
A 

effective activation energy for 
chemical diffusivity in alloy A, 
J/mole 

A 

Q„ effective activation energy for 

s 

solution of elements of alloy A, 
J/mole 
A 

Q effective activation energy for 

bl 

solution rate of elements of 
alloy A, J/mole 

0^ activation energy for viscous flow, 

J/mole 

q heat loss, J/hr 

R gas constant, J/(°K)(mole) 

o 

R m coefficient of multiple determina- 

tion 

r mass flux resistance, hr/m 

r Q radius of diffusing species 

S material-transport step 

s standard deviation 

s^ variance 

T absolute temperature, °K 

t time, hr 

v liquid flow velocity, m/hr 

v average film flow velocity, m/hr 

w mass transport, kg 
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constants 

d 

diffusional process 

e >V x j 


i 

transport step 

K 

thermal conductivity, 

k 

transport step parallel path 


J/(m)(hr)(°K) 

r 

derived from corrosion data 

V 

2 

kinematic viscosity, m /hr 

s 

rate-controlling transport 

1 

convective diffusion parameter 


step 

P 

O 

density, kg/m 

w 

wall 

Sub scripts: 




b 

bulk; boundary -region process 
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